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The  design  and  analysis  of  CPV  systems  require  knowledge  of  direct  normal  solar  irradiation  but 
ground-based  measurements  of  these  data  are  only  available  for  very  few  locations.  Nowadays, 
meteorological  databases  that  estimate  direct  normal  irradiation  from  satellite  images  and  other  data 
sources  are  used.  However,  values  provided  by  the  different  existing  databases  show  large  dispersion 
due  to  different  estimation  methods,  input  data  and  base  years.  In  this  paper,  a  simplified  method  for 
calculating  direct  normal  irradiation  is  presented.  It  has  been  obtained  from  previous  models  proposed 
by  several  authors.  One  of  its  advantages  is  that  it  only  requires  latitude  and  global  horizontal 
irradiation  as  input  data.  As  global  irradiation  is  easy  to  find  or  measure,  the  procedure  becomes  a 
useful  tool  in  renewable  energy  applications.  The  accuracy  of  this  method  is  similar  to  that  of  the 
existing  databases  and  it  is  able  to  easily  generate  a  mass  of  direct  normal  irradiation  data  for  different 
areas  worldwide. 
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1.  Introduction 

Photovoltaic  devices  are  at  a  mature  technological  stage  with  a 
long  and  wide  field  experience,  but  if  this  energy  source  wants  to 
compete  against  other  renewable  energy  sources  or  even  to  get 
into  the  traditional  energy  generation  system,  further  scientific 
progress  in  the  behavior  of  these  kinds  of  systems  is  necessary. 
Concentrator  Photovoltaic  (CPV)  technology  can  have  a  main  role 


Abbreviations:  CPV,  concentrator  photovoltaic;  DNI,  direct  normal  irradiance  or 
irradiation;  TMY,  typical  meteorological  year;  HC,  HelioClim;  HCPV,  high  con¬ 
centration  photovoltaic 
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in  this  challenge,  as  it  has  proved,  in  recent  published  research,  to 
have  the  potential  to  achieve  high  levels  of  energy  conversion 
performance  [1], 

CPV  is  based  on  the  use  of  optical  devices  that  increase  the 
light  received  on  the  solar  cell  surface.  The  concentration 
mechanism  is  achieved  by  lenses  or  mirrors  that  either  refract 
or  reflect  the  incoming  light  from  the  sun  on  top  of  the  cell 
receiver  surface.  Depending  on  the  concentration  process,  as  well 
as  on  the  concentration  area,  there  is  a  wide  range  of  optical 
device  configurations  that  can  be  used  in  the  implementation  of  a 
CPV  module  [2], 

The  use  of  concentrators  implies  that  CPV  systems  only  work 
with  the  Direct  Normal  Irradiance  (DNI).  So  it  is  necessary  to 
know  DNI  data  in  order  to  estimate  the  energy  that  will  be 
produced  by  the  system,  perform  economic  analysis,  supervise 
plant  operation,  etc.  However,  DNI  ground-based  measurements 
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Nomenclature 

B0  solar  constant  (W/m2) 

Sdn  solar  declination  of  the  day  in  the  middle  of  each 
month  (radians) 

dn  day  number  of  the  year  (dimensionless) 

eb  dm  eccentricity  correction  factor  of  the  earth's  orbit  of 
the  day  in  the  middle  of  each  month  (dimensionless) 

Gc  hm( 0)  monthly  average  hourly  global  irradiance  on  a 
horizontal  surface  (W/m2) 

GB  hm( 0)  monthly  average  hourly  direct  irradiance  on  a 
horizontal  surface  (W/m2) 

Gd  hm{ 0)  monthly  average  hourly  diffuse  irradiance  on  a  hor¬ 
izontal  surface  (W/m2) 

Gg  hm(S)  monthly  average  hourly  direct  normal  irradiance 
(W/m2) 


Hc  dm( 0)  monthly  average  daily  global  irradiation  on  a  hori¬ 
zontal  surface  (Wh/m2) 

HBO  dm( 0)  monthly  average  daily  extraterrestrial  irradiation  on 
a  horizontal  surface  (Wh/m2) 

HB  dm( 0)  monthly  average  daily  direct  irradiation  on  a 
horizontal  surface  (Wh/m2) 

HD  dm( 0)  monthly  average  daily  diffuse  irradiation  on  a 
horizontal  surface  (Wh/m2) 

HB  dm(S )  monthly  average  daily  direct  normal  irradiation 
(Wh/m2) 

Hb  a(S)  annual  direct  normal  irradiation  (Wh/m2) 

Kjm  clearness  index  (dimensionless) 

I<Dm  diffuse  fraction  of  global  radiation  (dimensionless) 

cp  latitude  (°) 

ca  solar  hour  angle  (radians) 

ms  dm  Sunrise  hour  angle  of  the  day  in  the  middle  of  each 
month  (radians) 


are  expensive  and  rarely  available  due  to  the  cost  and  sophistica¬ 
tion  of  measurement  devices  and  data  processing  requirements. 
There  is  a  particular  lack  of  data  in  the  Sunbelt  countries,  which 
are  more  favorable  for  the  use  of  CPV. 

This  lack  of  DNI  data  contrasts  with  the  availability  of  global 
horizontal  irradiation  data,  which  are  easy  to  find  or  measure.  So 
it  seems  to  be  interesting  to  have  a  procedure  that  obtains  DNI 
from  global  irradiation. 

We  can  find  in  the  literature  several  models  for  calculating  DNI 
from  different  parameters  [4-8].  These  models  are  accurate  but 
some  of  them  are  complex  and  difficult  to  apply,  especially 
because  they  require  input  data  that  is  not  easy  to  obtain,  for 
instance  measurements  taken  from  satellite  images.  In  this  paper, 
a  simplified  method  obtained  from  previous  models  is  presented. 
Its  aim  is  not  to  improve  the  existing  models,  but  to  easily  provide 
DNI  only  by  using  latitude  and  global  horizontal  irradiation  as 
input  data.  The  accuracy  of  the  method  seems  to  be  acceptable  for 
its  use  in  CPV  applications. 


2.  DNI  data  sources 

Preferably,  Typical  Meteorological  Year  (TMY)  data  should  be 
used  for  the  analysis  and  design  of  CPV  systems.  These  datasets 
are  based  on  measurements  taken  from  ground  meteorological 
stations  for  ten  years  [3],  However,  due  to  equipment  costs,  such 
measurements  are  scarce  and  TMY  datasets  are  available  for  only 
a  very  few  locations.  An  alternative  solution  to  this  problem  is  the 
use  of  a  model  that  derives  DNI  from  satellite  data.  In  recent 
years,  several  methods  for  estimating  direct  normal  irradiance 
and  irradiation  from  satellite  data  have  been  proposed  [4-8]. 
Nevertheless,  the  use  of  these  procedures  is  often  difficult  as  is 
obtaining  the  satellite  data.  Thus,  spatial  databases  that  directly 
provide  DNI  are  usually  utilized. 

There  are  a  number  of  spatial  databases  that  provide  DNI 
values  for  different  places  and  time  intervals  (Meteonorm  [9], 
NASA-SSE  [10],  PVGIS  [11],  PVSAT  [12],  Satel-Light  [13],  SoDa 
[14],  SOLEMI  [15],  etc.).  These  databases  use  different  kinds  of 
input  data  (satellite  data,  ground-based  measurements)  and 
different  procedures  in  order  to  estimate  DNI  values.  They  are  a 
very  useful  tool  for  the  design  and  evaluation  of  CPV  systems. 
However,  some  databases  show  weaknesses:  there  is  a  lack  of 
continuity  in  the  data  because  the  information  comes  from  a 
limited  number  of  ground  stations  and  there  are  no  hourly  DNI 
data,  except  for  a  very  few  sites.  Although  the  main  problem  is 


uncertainty  because  DNI  values  provided  by  the  different  data¬ 
bases  show  large  dispersion.  The  result  is  a  fragmentation  of 
services,  each  having  its  own  mechanism  of  access  and  all  giving 
different  outputs  due  to  different  methods,  input  data  and 
base  years. 

In  order  to  analyze  this  uncertainty,  a  comparison  between 
measured  DNI  values  taken  from  ground  meteorological  stations 
located  at  three  Spanish  places  (Madrid,  40.41 7N  -3.700W; 
Murcia,  37.983N  1.130W  and  Jaen,  37.766N-3.790W)  and  DNI 
values  provided  by  three  databases  (SoDa,  NASA-SSE  and  PVGIS) 
for  these  three  locations  has  been  made. 

DNI  measurements  for  Madrid  have  been  carried  out  by  the 
Solar  Energy  Institute  (Instituto  de  Energia  Solar,  IES)  in  2006, 
2007  and  2010.  DNI  measurements  for  Murcia  have  been  carried 
out  by  the  Meteorology  State  Agency  (Agencia  Estatal  de  Meteor- 
ologia,  AEMET)  from  2005  to  2010.  And  DNI  measurements  for 
Jaen  have  been  carried  out  by  the  University  of  Jaen  in  2001, 2002, 
2003  and  2010.  In  every  case,  a  pyrheliometer  with  wavelength 
range  from  0.2  to  4  pm  has  been  used  for  the  DNI  measurements. 

SoDa  website  disseminates  the  HelioClim  (HC)  databases 
which  contain  radiation  values  at  ground  level  obtained  from 
the  processing  of  the  images  taken  by  Meteosat  satellites.  In  this 
case,  the  HC  databases  are  processed  by  MINES  ParisTech.  On  the 
other  hand,  NASA-SSE  release  6  combines  results  from  GEWEX/ 
SRB  3  and  ISCCP  projects  with  NCAR  reanalysis  products  (1983- 
2005).  The  database  directly  provides  monthly  and  annual  DNI 
values  around  the  world.  Finally,  PVGIS  includes  solar  radiation 
developed  by  a  combination  of  solar  radiation  models  and 
interpolated  ground  observations,  representing  the  period 
1981-1990.  It  does  not  directly  provide  DNI  but  this  can  be 
estimated  from  global  and  diffuse  radiation  components. 

Table  1  compares  measured  values  taken  from  the  meteor¬ 
ological  stations  to  values  provided  by  the  databases.  Annual  DNI 
value  measured  at  Madrid  is  1882  kWh/m2  while  the  different 
databases  give  values  between  1814  and  1964  kWh/m2.  This 
implies  an  error  between  —3.6%  and  4.4%.  Annual  DNI  value 


Table  1 

Annual  DNI  values  for  three  Spanish  locations  (kWh/m2). 


Location 

Measured 

SoDa 

NASA-SSE 

PVGIS 

MADRID 

1882 

1947 

1964 

1814 

MURCIA 

1961 

1738 

1905 

1870 

JAEN 

1984 

2208 

2132 

1745 
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measured  at  Murcia  is  1961  l<Wh/m2  while  the  different  data¬ 
bases  give  values  between  1738  and  1905  kW  h/m2.  This  implies 
an  error  between  -11.4%  and  -2.9%.  And  annual  DNI  value 
measured  at  Jaen  is  1984  kW  h/m2  while  the  different  databases 
give  values  between  1745  and  2208  l<Wh/m2.  This  implies  an 
error  between  -12.0%  and  11.3%.  Fig.  1  shows  the  annual  DNI 
error  for  the  three  databases.  A  clear  trend  cannot  be  observed, 
although  PVGIS  underestimates  DNI  for  these  three  locations  and 
the  three  databases  underestimate  Murcia's  DNI  value.  As  a 
conclusion,  it  can  be  assumed  that  it  is  difficult  to  get  reliable 
DNI  data  because  each  database  gives  different  values  and  thus, 
there  are  sources  of  uncertainty. 

Other  authors  have  reported  the  same  conclusions  about  uncer¬ 
tainty  [16,17],  Suri  et  al.  [16]  cross-compared  maps  of  annual  DNI 
from  five  databases  offering  solar  resource  and  climate  data  for 
Europe:  Meteonorm  [18],  NASA-SSE  [19],  PVGIS  [20],  Satel-Light  [21] 
and  SOLEMI  [22].  From  these  databases,  the  following  maps  were 
calculated:  (1)  overall  average  estimated  by  a  simple  averaging, 
(2)  five  maps  of  differences  to  the  average,  and  (3)  standard  devia¬ 
tion.  While  maps  of  differences  show  the  deviation  of  the  particular 
dataset  from  the  overall  average,  the  standard  deviation  indicates 
user's  uncertainty  given  by  the  magnitude  of  differences  between  the 
combined  data  sources.  Results  show  that  90%  of  the  analyzed  area 
presents  an  uncertainty  in  annual  DNI,  expressed  by  relative  stan¬ 
dard  deviation,  of  up  to  17%.  This  high  uncertainty  indicates  that  it  is 
necessary  to  improve  the  databases  and  the  radiation  models.  Some 
efforts  have  been  made  in  this  field.  For  instance,  the  MESoR  project 
[17],  funded  by  the  European  Commission,  reduces  the  associated 
uncertainty  by  setting  up  standard  benchmarking  rules  and  mea¬ 
sures  for  comparing  the  databases,  user  guidance  to  the  application 
of  resource  data  and  unifying  access  to  various  databases. 

Below,  a  very  simple  procedure  for  the  calculation  of  DNI 
based  on  previous  existing  models  is  presented.  It  can  be  easily 
implemented  in  a  spreadsheet  or  in  computer  applications  in 
renewable  energy.  The  aim  of  the  proposed  procedure  is  not  to 
provide  high  accuracy.  It  tries  to  meet  the  basic  requirements  for 
the  analysis  of  CPV  systems  by  providing  DNI  errors  similar  to 
those  of  the  different  meteorological  databases.  Inputs  of  the 
procedure  are  the  global  horizontal  irradiation  values  which, 
unlike  DNI  values,  are  often  measured  in  the  meteorological 
stations  and  are  easy  to  find.  The  procedure  does  not  improve 
the  existing  methods  but  is  an  easy-to-use  choice  useful  for 
teaching  and  providing  acceptable  accuracy. 


In  what  follows,  according  to  [23,24],  letter  G  will  be  used  to 
indicate  irradiance  or  the  power  incident  on  a  surface  per  unit 
area  of  surface  (units:  W/m2)  and  letter  H  will  be  used  to  indicate 
irradiation  or  the  energy  incident  on  a  surface  over  a  specified 
period  of  time  per  unit  area  of  surface  (units:  W  h/m2). 

At  the  high  concentration  rates  involved  in  High  Concentration 
Photovoltaic  (HCPV)  systems,  time  resolution  as  well  as  uncer¬ 
tainty  in  radiation  data  play  a  complex  role  [25].  However,  some 
authors  have  proven  that  the  complexity  of  a  solar  radiation 
model  needed  for  yearly  energy  calculations  is  very  low  [26]. 
Twelve  values  of  monthly  mean  solar  irradiation  are  enough  to 
estimate  yearly  energy  with  small  errors.  Also,  time  resolutions 
shorter  than  hourly  samples  do  not  significantly  improve  the 
energy  estimation  results. 

So,  the  start-up  data  required  for  the  calculations  are  the  site 
latitude,  tp,  and  the  twelve  values  of  monthly  average  daily  global 
irradiation  on  a  horizontal  surface,  Hcdm( 0).  Results  of  the  method 
are  the  monthly  average  direct  normal  irradiance  at  one  hour 
intervals  for  the  twelve  months  of  the  year,  CB  hm(S),  the  twelve 
values  of  monthly  average  daily  direct  normal  irradiation,  HB  dm{S), 
and  the  value  of  the  annual  direct  normal  irradiation,  HB  a(S). 
The  proposed  method  performs  the  sequence  shown  in  Fig.  2  in 
order  to  calculate  the  DNI  values. 

The  method  is  based  on  easy-to-use  equations  proposed  by 
several  authors  [27]  which  provide  a  good  approximation  for  a  wide 
range  of  latitudes.  Other  more  accurate  equations  and  correlations 
[28-30]  have  not  been  taken  into  account  because  of  their  complex¬ 
ity  or  because  they  are  only  valid  for  a  specific  area. 

The  proposed  method  is  based  on  the  following  procedure: 


I.  Solar  declination  of  the  day  in  the  middle  of  each  month,  in 
radians,  5dn  [31]: 


= 


23.4571 

180 


sin 


d„  +  284\ 
365  J 


(1) 


where  dn  is  the  day  number  of  the  year,  ranging  from  1  on 
January  1  to  365  on  December  31.  The  values  corresponding 
to  January  15,  February  15,  March  15,  etc.  will  be  used. 

II.  Eccentricity  correction  factor  of  the  earth’s  orbit  of  the  day  in 
the  middle  of  each  month,  dimensionless,  e0  dm  [32]: 

£o  dm  =  1+0.033  x  cos  (^27t  x  (2) 


3.  Description  of  the  proposed  method 

In  this  section,  a  procedure  that  calculates  DNI  values  from 
global  horizontal  irradiation  values  is  presented. 


■  MADRID 

□  MURCIA 

□  JAEN 


III.  Sunrise  hour  angle  of  the  day  in  the  middle  of  each  month,  in 
radians,  cos  dm  [33]: 

cosdm  =  -cos_1(-tan  5dn  x  tan  tp)  (3) 


1.  Calculation  of  daily  extra  terrestrial  irradiation  ona  horizontal  surface. 

2.  Calculation  of  daily  diffuse  irradiation  on  a  horizontal  surface. 

3.  Diffuse  irradiance  calculation  from  irradiation. 

4.  Global  irradiance  calculation  from  irradiation. 

5.  Calculation  of  hourly  direct  normal  irradiance. 

6.  Calculation  of  daily  direct  normal  irradiation. 

7.  Calculation  of  annual  direct  normal  irradiation. 


Fig.  1.  Annual  DNI  error  for  the  selected  locations  and  databases. 


Fig.  2.  Flow  diagram  of  the  proposed  method. 
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IV.  Monthly  average  daily  extraterrestrial  irradiation  on  a  hor¬ 
izontal  surface  of  each  month,  in  W  h/m2,  HBO  dm(0)  [34]: 

24 

HB0  dm(G)  =  —B0  X  Eo  dm  X  (COS(/>  ■  COS(5dn) 

x(ttAdm  x  COS OJsdm  sinC(2sdm)  (4) 

being  60=1367  W/m2  the  solar  constant  [35].  It  is  assumed 
that  the  extraterrestrial  irradiation  of  the  day  in  the  middle 
of  each  month  almost  equals  the  monthly  average  daily 
extraterrestrial  irradiation  of  this  month. 

V.  Monthly  average  clearness  index,  dimensionless,  I<Tm  [36]: 

^Tm  =  Hcdm(O) /Hfiodm(O)  (5) 


XIV.  Annual  direct  normal  irradiation,  in  W  h/m2,  HB  a(S): 

12 

Hb  a(S)  ^J2dmiX  HB  dm  i(S)  (16) 

i  =  1 

where  dm,  is  the  number  of  days  for  the  “i"  month  and  HB  dm  ,(S) 
is  the  monthly  average  daily  direct  normal  irradiation  for  the 
“i”  month. 


4.  Results  and  discussion 


VI.  Monthly  average  diffuse  fraction  of  global  radiation,  dimen¬ 
sionless,  I<Dm  [37]: 

XDm  =  l-1.13XTm  (6) 

VII.  Monthly  average  daily  diffuse  irradiation  on  a  horizontal 
surface  of  each  month,  in  Wh/m2,  HD  dm( 0): 

Hd  dm(0)  —  ffcdm(O)  x  K Dm  (7) 

VIII.  Monthly  average  daily  direct  irradiation  on  a  horizontal 
surface  of  each  month,  in  W  h/m2,  HB  dm(0): 

Hb  dm(0)  =  Hg  dm(0)—HD  dm( 0)  (8) 


IX.  Monthly  average  hourly  diffuse  irradiance  on  a  horizontal 
surface,  in  W/m2,  CD  flm(0)  [38]: 


CD„m(0)  =  HDdm(0)^x  (- - C0S“C°SC°\d: ”  - 

24  \&)s  dm  x  COS COsdm  SinCOs 


(9) 


where  m  is  the  solar  hour  angle  in  radians  (-7t  for  0:00 
hours  and  +7t  for  24:00  hours).  Irradiance  values  are 
calculated  at  one  hour  intervals  (Aco=27t/24  rad). 

Monthly  average  hourly  global  irradiance  on  a  horizontal 
surface,  in  W/m2,  Cc  hm{ 0)  [39]: 


Gg  hm(0)  =  Hc  dm(0)^(a+b 


x  COSOJ) 


COSO)-COStUsdm 


(10) 


Wh  dm  x  COS (Dsdm  SintOs  dmj 
where  a  and  b  parameters  are  calculated  by  the  formulae: 
a  =  0.4090-0.5016  sin(ms  dm  + 1.047)  (11) 


b  =  0.6609  +  0.4767sin(cos 


-1.047) 


(12) 


In  this  section,  the  proposed  method  is  applied  to  three 
Spanish  locations  (Madrid,  Murcia  and  Jaen)  in  order  to  estimate 
both  the  monthly  average  daily  direct  normal  irradiation  of  each 
month,  Hb  dm(S ),  and  the  annual  direct  normal  irradiation,  HB  „(S). 
These  locations  correspond  to  the  locations  analyzed  in  Section  2. 
They  have  been  selected  because  ground-based  measurements  of 
annual  DNI  are  available  and  thus,  a  comparison  between  mea¬ 
sured  and  calculated  values  can  be  made.  In  what  follows, 
irradiation  values  are  expressed  in  kW  h/m2. 

Table  2  shows  the  input  data  used  for  the  calculations,  i.e„  the 
site  latitude,  cp,  and  the  twelve  values  of  monthly  average  daily 
global  irradiation  on  a  horizontal  surface,  Hc  dm( 0).  The  global 
irradiation  values  have  been  obtained  from  PVGIS. 

With  these  input  data,  the  procedure  described  in  Section  3  can 
be  easily  applied.  Results  are  shown  in  Table  3,  including  the  twelve 
values  of  monthly  average  daily  direct  normal  irradiation,  HB  dm(S), 
and  the  value  of  the  annual  direct  normal  irradiation,  HB  a(S). 

In  order  to  compare  measured  and  calculated  DNI  values,  the 
annual  DNI  error  has  been  computed  for  each  location: 

error(%)  =  Hb  Wufarerf-H8  a(S )measured  x  1 Q0  (17) 

"B  aWJmeasured 

Table  4  shows  the  results.  The  proposed  method  overestimates 
DNI  for  the  three  locations.  Madrid’s  error  is  3.9%  while  the  three 
meteorological  databases  selected  in  Section  2  (SoDa,  NASA-SSE 
and  PVGIS)  gave  errors  between  -3.6%  and  4.4%.  Murcia's  error  is 
1.2%  while  the  databases  errors  were  between  -11.4%  and 
-2.9%.  And  Jaen's  error  is  0.2%  while  the  databases  errors  were 
between  -12.0%  and  11.3%.  As  a  conclusion,  annual  DNI  values 
produced  by  the  proposed  method  are  comparable  to  those  of  the 
existing  meteorological  databases.  Thus,  an  alternative  easy-to- 
use  method  can  be  used  to  compute  DNI. 


XI.  Monthly  average  hourly  direct  irradiance  on  a  horizontal 
surface,  in  W/m2,  CB  hm( 0): 

Gb  hm(0)  =  Gc  hm(0)— GD  /,m(0)  (13) 


XII.  Monthly  average  hourly  direct  normal  irradiance,  in  W/m2, 
Gb  US)  [40]: 


GBhm(S)  = 


GB  hm( 0) 


sin<5d„  x  sin<p  +  coS(5dn  x  cosi p  x  cosco 


(14) 


XIII.  Monthly  average  daily  direct  normal  irradiation,  in  Wh/m2, 
HB  dm(S): 

24 

HB  dm(S)  ^  ^2  GB  hmj(S )  x  (15) 

j=1 


where  CB  hmj(S)  is  the  monthly  average  hourly  direct  normal 
irradiance  at  “j”  hour  and  Af=l  h. 


Table  2 

Input  data  required  by  the  proposed  method  for  the  three  selected  locations: 
latitude  and  monthly  average  daily  global  irradiation  on  a  horizontal  surface. 


Location: 

Madrid 

Murcia 

Jaen 

latitude,  °: 

40.417 

37.983 

37.766 

Month 

He  dm(0).  kW  h/m2 

j 

1.98 

2.44 

2.46 

F 

2.69 

3.21 

3.16 

M 

4.44 

4.54 

4.65 

A 

5.09 

5.40 

5.29 

M 

6.50 

6.48 

6.59 

j 

7.23 

7.00 

7.19 

j 

7.33 

7.05 

7.07 

A 

6.44 

6.15 

6.22 

S 

4.98 

5.00 

4.94 

O 

3.36 

3.74 

3.74 

N 

2.14 

2.47 

2.49 

D 

1.61 

2.16 

2.09 

P.J.  Perez-Higueras  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  16  (2012)  5529-5534 


5533 


5.  Application 

The  simplified  procedure  described  in  this  paper  can  easily  be 
used  to  obtain  DNI  data  worldwide.  The  only  requirement  is  to 
know  global  horizontal  irradiation  values  for  the  area  of  interest. 
As  these  values  are  easy  to  find  or  measure,  the  proposed  method 


Table  3 

Output  data  obtained  by  the  proposed  method:  monthly  average  daily  DNI  of  each 
month  and  annual  DNI. 


Location: 

Madrid 

Murcia 

Jaen 

Month 

Hb  dm( s),  l<w  h/m2 

j 

3.12 

4.02 

4.02 

F 

3.44 

4.29 

4.10 

M 

5.72 

5.44 

5.69 

A 

5.04 

5.39 

5.14 

M 

6.75 

6.55 

6.78 

j 

7.86 

7.19 

7.60 

j 

8.37 

7.61 

7.65 

A 

7.55 

6.58 

6.72 

S 

6.20 

5.79 

5.62 

0 

4.56 

5.03 

4.98 

N 

3.21 

3.64 

3.65 

D 

2.44 

3.73 

3.42 

a(S),  kW  h/m2 

Annual 

1955 

1985 

1988 

Table  4 

Comparison  between  measured  annual  DNI  and  annual  DNI  calculated  by  the 
proposed  method. 


Location  HB  a(S),  kW  h/m2 


Measured 

Calculated 

Error  (%) 

Madrid 

1882 

1955 

3.9 

Murcia 

1961 

1985 

1.2 

Jaen 

1984 

1988 

0.2 

becomes  a  very  useful  tool  in  renewable  energy  applications, 
particularly  in  CPV  applications. 

As  an  example  of  such  applications,  a  complete  map  of  the 
Spanish  annual  DNI  is  presented  in  Fig.  3.  It  has  been  obtained 
from  global  irradiation  data  of  the  47  peninsular  provincial 
capitals  of  Spain.  For  each  location,  annual  DNI  has  been  calcu¬ 
lated  by  using  a  Microsoft  Excel  spreadsheet  with  the  procedure 
described.  This  spreadsheet  is  available  at  http://www.ujaen.es/ 
investiga/solar/dni.xlsx.  The  DNI  data  obtained  has  been  graphi¬ 
cally  represented  with  the  help  of  the  Surfer  program  by  Golden 
Software  [41]. 

This  example  proves  that  the  proposed  method  can  be  used  to 
obtain  a  mass  of  DNI  data  for  different  regions  worldwide  with 
acceptable  accuracy. 


6.  Conclusions 

DNI  data  provided  by  the  existing  meteorological  databases  is 
obtained  from  satellite  images  and  ground-based  measurements 
by  using  different  methods,  input  data  and  base  years.  While 
these  databases  are  a  useful  tool  in  renewable  energy  applica¬ 
tions,  they  present  uncertainty. 

A  simplified  procedure  for  estimating  DNI  from  global  irradia¬ 
tion  based  on  previous  existing  models  has  been  developed.  It 
provides  accuracy  similar  to  that  of  the  existing  databases.  As 
global  irradiation  data  are  easy  to  find  or  measure,  the  method 
becomes  a  good  choice  for  the  design  and  analysis  of  CPV  and 
other  renewable  energy  systems.  It  can  be  easily  implemented  in 
computer  programs  and  is  useful  for  teaching. 

The  accuracy  of  the  method  has  been  shown  by  comparing 
calculated  and  measured  DNI  for  three  Spanish  locations.  Also,  the 
estimation  error  has  been  compared  with  the  errors  provided  by 
three  existing  databases.  Results  indicate  that  the  procedure  can 
be  used  as  an  alternative  to  other  sources  of  DNI  data. 

In  order  to  prove  the  ability  of  the  proposed  method  to  easily 
obtain  a  mass  of  DNI  data  for  different  areas  worldwide,  a  map  of 
the  Spanish  annual  DNI  has  been  developed  and  presented. 


(C)  Universidad  de  Jaen 
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Fig.  3.  Map  of  the  Spanish  annual  DNI  provided  by  the  proposed  method. 
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